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o 
o 

i— i Quantum computers have the potential to solve certain interesting problems significantly faster than 

T—H classical computers. To exploit the power of a quantum computation it is necessary to perform inter- 

00 qubit operations and generate entangled states. Spin qubits are a promising candidate for implement- 

CM 

00 m g a quantum processor due to their potential for scalability and miniaturization. However, their 

weak interactions with the environment, which leads to their long coherence times, makes inter- qubit 
operations challenging. We perform a controlled two-qubit operation between singlet- triplet qubits 
using a dynamically decoupled sequence that maintains the two-qubit coupling while decoupling each 

> 

. !pH qubit from its fluctuating environment. Using state tomography we measure the full density matrix of 

the system and determine the concurrence and the fidelity of the generated state, providing proof of 
entanglement. 

Singlet-triplet (S-Tb) qubits, a particular realization of spin qubits [1-7], store quantum information in the joint 
spin state of two electrons [8-10]. The basis states for the S-T qubit can be constructed from the eigenstates 
of a single electron spin, | f> and | |). We choose \S) = ± (| U> - | |t» and \T ) = -±= (I t i) + I it)) for the logical 
subspace of the S- To qubit because these states are insensitive to uniform fluctuations in the magnetic field. 
The qubit can then be described as a two level system with a representation on a Bloch sphere shown in Fig. la. 
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Universal quantum control is achieved using two physically distinct operations that drive rotations around the 
x and z-axes of the Bloch sphere [11]. Rotations around the z-axis of the Bloch sphere are driven by the exchange 
splitting, /, between \S) and | To), and rotations around the x-axis are driven by a magnetic field gradient, /\B Z 
between the electrons. 

We implement the S-T qubit by confining two electrons to a double quantum dot (QD) in a two 
dimensional electron gas (2DEG) located 91nm below the surface of a GaAs-AlGaAs heterostructure. We deposit 
local top gates using standard electron beam lithography techniques in order to locally deplete the 2DEG and 
form the QDs. We operate between the states (0,2) and (1,1) where {ni y nR) describes the state with n^TiR) 
electrons in the left (right) QD. The |S) and \T ) states, the logical subspace for the qubit, are isolated by 
applying an external magnetic field of B =700mT in the plane of the device such that the Zeeman splitting makes 
T+ = | 1 1), and T- = | j j) energetically inaccessible. The exchange splitting, /, is a function of the difference in 
energy, e, between the levels of the left and right QDs. Pulsed DC electric fields rapidly change e, allowing us to 
switch / on, which drives rotations around the z-axis. When / is off the qubit precesses around the x-axis due to 
a fixed A5 Z , which is stabilized to ^zlin =30MHz by operating the qubit as a feedback loop between interations 
of the experiment[12] . Dephasing of the qubit rotations reflects fluctuations in the magnitude of the two control 
axes, / and A5, caused by electrical noise and variation in the magnetic field gradient, respectively. The qubit 
is rapidly (<50ns) initialized in |S) by exchanging an electron with the nearby Fermi sea of the leads of the QD 
in a region of (0,2) where only |S) is accessible, and the qubit state is read out using standard Pauli blockade 
techniques, where e is quickly tuned to the regime where S occupies (0,2) and T occupies (1,1), allowing the 
qubit state to be determined by the proximal charge sensor. The charge state of the qubit is rapidly read (~ Ifis) 
using standard RF-techniques [13, 14] on an adjacent sensing QD. 

In order to make use of the power of quantum information processing it is necessary to perform two qubit 
operations in which the state of one qubit is conditioned on the state of the other [15]. To investigate two- 
qubit operations we fabricate two adjacent S-To qubits such that they are capacitively coupled, but tunneling 
between them is suppressed (Fig. lb). A charge sensing QD next to each qubit allows for simultaneous and 
independent projective measurement of each qubit (supplement). We use the electrostatic coupling between 
the qubits to generate the two-qubit operation [16]. When / is nonzero, the S and T states have different charge 
configurations in the two QDs due to the Pauli exclusion principle (Fig. lc). This charge difference, which is a 
function of e, causes the |S) and | To) states in one qubit to impose different electric fields on the other qubit. 
Since / is a function of the electric field, the change imposed by the first qubit causes a shift in the precession 
frequency of the second qubit. In this way the state of the second qubit may be conditioned on the state of the 
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first qubit. More precisely, when a single qubit evolves under exchange, there exists a state-dependent dipole 
moment, d> between \S) and |To) resulting from their difference in charge occupation of the QDs. Therefore, 
when simultaneously evolving both qubits under exchange, they experience a capacitively mediated, dipole- 
dipole coupling that can generate an entangled state. The two-qubit Hamiltonian is therefore given by: 



h I hi 

^2-qubi t=-\MVz®D+MI® Oz) + — {iP Z -D*iPz-l))+ (0" x /) + ^B Z} 2 (1 ® O x ) 



(1) 



where o X} y }Z are the Pauli matrices, J is the identity operator, AB z j and J\ are the magnetic field gradients and 
the exchange splittings (i=l,2 for the two qubits), and 7i2 is the two-qubit coupling, which is proportional to the 
product of the dipole moments in each qubit. For a two level system with constant tunnel coupling, the dipole 
moment scales as d\ oc Empirically, we find that for experimentally relevant values of //, ^ oc ]\ (e), so that 
hi oc /1/2. As with the single qubit operations, this two-qubit operation requires only pulsed DC electric fields. 

In principle, evolving both qubits under exchange produces an entangling gate. However, the time 
to produce this maximally entangled state exceeds the inhomogeneously broadened coherence times of 
each individual qubit, rendering this simple implementation of the two-qubit gate ineffective. To mitigate 
this we use a dynamically decoupled entangling sequence[17, 18] (Fig. Id). In this sequence, each qubit 
is prepared in |S) and is then rotated by | around the x-axis Ui =0, /\B Zf il2n ^30MHz) to prepare a 
state in the x-y plane. The two qubits are subsequently both evolved under a large exchange splitting 
(h/2n -280MHZ, J 2 /2n « 320MHz » A£ z ) for a time |, during which the qubits begin to entangle and 
disentangle. A 7T-pulse around the x-axis [AB Z ) is then applied simultaneously to both qubits, after which the 
qubits are again allowed to exchange for a time |. This Hahn echo-like sequence[19] removes the dephasing 
effect of noise that is low frequency compared to 1/r, while the 7i-pulses preserve the sign of the two-qubit 
interaction. The resulting operation produces a CPHASE gate, which, in a basis of {|SS>, \TqS), \ST ), |r r )}, 
is an operation described by a matrix with e~ l6/2 ,l,l f e~ lG/2 on the diagonals. For r = r ent = ^ (0 = n) the 
resulting state is a maximally entangled generalized Bell state Y^ent) = g z;r(/<8)a > /+a ^ <8)/ ) /8 |^_), which differs from 
the Bell state |Y_> = ± (|SS) - | T T )) by single qubit rotations. 

In order to characterize our two-qubit gate and verify that we produce an entangled state we perform 
two-qubit state tomography and extract the density matrix and appropriate entanglement measures. The 
tomographic procedure is carefully calibrated with minimal assumptions in order to avoid adding spurious 
correlations to the data that may artificially increase the measured degree of entanglement (supplement). We 
choose the Pauli set representation of the density matrix[15, 20, 21], where we measure and plot the 16 two- 
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qubit correlators (ij) = (OiOj) where <7; are the Pauli matrices and i,j e {I,X,Y,Z}. As a first measure of 
entanglement, we evaluate the concurrence [22] (Fig. 2a), C (p) = max{0, X4-X3-X2- X\} for different r, where 
p is the experimentally measured density matrix, and X\ are the eigenvalues, sorted from largest to smallest, 
of the matrix R = yj^fppy/py and p = {a y <8> a y )p*{a y <8> a y ), and p* is the complex conjugate of p. A positive 
value of the concurrence is a necessary and sufficient condition for demonstration of entanglement [22]. For 
r = 140ns we extract a maximum concurrence of 0.44. 

While a positive value of the concurrence is a definitive proof of entanglement, it alone does not verify that 
the two-qubit operation produces the intended entangled state. In order to better characterize the generated 
quantum state, we evaluate another measure of entanglement, the Bell state fidelity, F = ' en t\pV¥ ent) • This 
may be interpreted as the probability of measuring our two-qubit state in desired Went)- Additionally, for all 
non-entangled states one can show that F < 0.5 [23, 24]. In terms of the Pauli basis, the Bell state fidelity takes 
the simple form F = \P ent • F 'experiment where P ent and P experiment are the Pauli sets of a pure target Bell state 
and of the experimentally measured state, respectively. For our target state VV e nt)y the resulting Pauli set is given 
by (XZ) = (ZX) = (YY) =1, with all other elements equal to zero (Fig 3a). 

In an idealized, dephasing-free version of the experiment, as r increases and the qubits entangle and 
disentangle, we expect the nonzero elements of the Pauli set for the resulting state to be 

(YI) = (IY) = cos(/i 2 t), (XZ) = (ZX) = sin(/i 2 T), (YY) = 1 (2) 

Dephasing due to electrical noise causes the amplitudes of the Pauli set to decay. However, the two-qubit 
Hamiltonian (equation(l)) includes rapid single-qubit rotations around the S-Tb axis U1J2 y>]\2l2n « 1MHz) 
that change with r due to imperfect pulse rise times in the experiment. These add additional single-qubit 
rotations around the S-T axis of each qubit, which are not accounted for in equation (2). We determine the 
angle of the single qubit rotations by performing a least-squares fit of the experimental data to modified form of 
equation (2) that accounts for these rotations and dephasing. The decays due to dephasing are fit by calculating 
p (f) in the presence of noise on 7i and /2, which leads to decay of certain terms in the density matrix [25, 26]. 
For the present case where J12 « J\> h> we neglect the two-qubit dephasing, which is smaller than single-qubit 
dephasings by a factor of j^, -300, and we extract a separate dephasing time for each individual qubit. We 
remove the single-qubit rotations numerically in order to simplify the presentation of the data (Fig. 3e). The 
extracted angles exhibit a smooth monotonic behavior which is consistent with their underlying origin (see 
supplemental information). 
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In the absence of dephasing we would expect the Bell state fidelity to oscillate between 0.5 for an 
unentangled state and 1 for an entangled state as a function of r. This oscillation is caused by the phase 
accumulated by a CPHASE gate between the two qubits. However, the qubits dephase as the state becomes 
increasingly mixed, and this oscillation decays to 0.25. Indeed, this behavior is observed (Fig. 2b) : for very short 
r there is very little dephasing present, and the qubits are not entangled. As r increases the Bell state fidelity 
increases as the qubits entangle, reaching a maximum value of 0.72 at r = 140ns. As r is increased further, we 
continue to see oscillations in the Bell state fidelity , though due to dephasing, they do not again rise above 0.5. 

Fig. 2c shows these oscillations in Bell state fidelity as a function of r for several different values of / as 
e is changed symmetrically in the two qubits. We see that as the value of / increases in the two qubits, the 
time required to produce a maximally entangled state, r ent , decreases, but the maximum attainable fidelity is 
approximately constant. This is consistent with the theory that Juoc ^ • ^ oc 7i • J2. 

To further understand the evolution of the quantum state, we focus on one value of / and compare the 
measured Pauli set to that expected from single-qubit dephasing rates and J12 (see supplemental information). 
Fig. 3a shows the Pauli set for the measured and expected quantum states for r =40ns, which shows three large 
bars in the (YI),(IY), and (YY) components of the Pauli set. This is a nearly unentangled state. At t =140ns, 
we see weight in the in the (XZ), (ZX), (YY) components of the Pauli set (Fig. 3b), and we extract a Bell state 
fidelity of 0.72, which demonstrates the production of an entangled state. For r = r ent = ^= 160ns (Fig. 3c) we 
see a similar state to r = 140ns, but with less weight in the single qubit components of the Pauli set. This state 
corresponds to the intended CPHASE of n, though it has a slightly lower fidelity than the state at r = 140ns due 
to additional decoherence. Finally, at t = =320ns (Fig. 3d), where we expect the state to be unentangled, 
we again see large weight in the (YI), (IY), and (YY) components of the Pauli set, though the bars are shorter 
than the Pauli set for r =40ns, due to dephasing of the qubits. We plot the entire Pauli set as a function of time 
(Fig. 3e), which clearly shows the predicted oscillation (equation(2)) between (YI) y (IY) and (XZ), (ZX), with 
decays due to decoherence. 

The two-qubit gate that we have demonstrated is an important step toward establishing a scalable 
architecture for quantum information processing in S-T qubits. State fidelity is lost to dephasing from 
electrical noise, and decreasing the ratio ^g^, where T^ cho is the single-qubit coherence time with an echo 
pulse, is therefore paramount to generating high-fidelity Bell states. Large improvements can be made by 
introducing an electrostatic coupler between the two qubits [27] in order to increase the two-qubit coupling 
(/12) and reduce r ent . We estimate that in the absence of other losses, if an electrostatic coupler were 
used, a Bell state with fidelity exceeding 90% could be produced. Additional improvements can be made 
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by studying and mitigating the origins of charge noise to increase T£ c . This would allow future tests of 
complex quantum operations including quantum algorithms and quantum error correction. The generation 
of entangled also states opens the possibility of studying the complex dynamics of the nuclear environment, 
which is a fundamental, quantum, many body problem. 
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a |S> b 




Figure 1 1 Two-qubit coupling scheme, a, A Bloch sphere can be used to describe the states of the effective two-level system defined by 
the singlet and triplet states of the qubit, with the z-axis along the S- To axis and the x-axis along the 1 1 J)/| J |> axes, b, An SEM image of 
the top of the device used shows gates used to define the S-Tq qubits (white), dedicated ns control leads, the approximate locations of 
the electrons in the two qubits (red), and current paths for the sensing dots (green arrows). The left qubit uses the LR and LL electrons, 
while the right qubit uses the RL and RR electrons, c, A schematic of the electronic charge configurations for the |S> (blue) and the | Tq) 
(red) at non-zero /. This difference in charge configuration is the basis for the electrostatic coupling between the qubits d, The pulse 
sequence used to entangle the qubits: initialize each qubit in a |S>, perform a nl2 rotation around the x-axis, allow the qubits to evolve 
under exchange for a time t/2, perform a 7r-rotation around the x-axis, thereby decoupling the qubits from the environment but not 
each- other, evolve under exchange for t/2, and perform state tomography to determine the resulting density matrix (see supplemental 
information.) 
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Figure 2 | Proof of entanglement: concurrence and state fidelity, a, A plot of the difference of the sorted eigenvalues of the matrix R, 
which for positive values is equal to the concurrence c (p). States with a concurrence greater than zero (shaded region) are necessarily 

entangled, b, The fidelity with which the measured state approximates the target Went) (blue), and e in ( (T y® I+I ® (T y} /4 \ x ¥ en t) (green), 
which differs from Went) by single qubit rotations and is the expected state for t = 37r/2/ 12 . The fidelity with which the measured state 
approximates a dephasing-free model of the entangling operation (red) shows smooth decay due to decoherence. The solid lines are 
fits to the data. Inset: The time to produce a maximally entangled state as a function of the change in e (and therefore /) in the two 
qubits. As / increases r e nt decreases (red), but the maximum attainable fidelity (green) is approximately constant. Arrows indicate 
which y-axis is to be used, c, The Bell state fidelity as a function of time for different values of / (offset) with guides to show where the 
fidelity exceeds 1/2 for each curve. As / increases in the two qubits the time to produce an entangled state, r e nt> decreases. 
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Figure 3 | Pauli set representation, a, The elements of the Pauli set for the measured density matrix and the state expected from the 
entangling and dephasing rates for small r (unentangled). b, The Pauli set of the measured and expected states for r = 140ns, which 
produces a maximum Bell state fidelity of 0.72. c, The Pauli set of the measured and expected states for t = r en t = 160ns, which is a 
CPHASE of 7i but does not the highest fidelity due to dephasing. d, The Pauli set for the measured and expected states for t =320ns, 
which is an unentangled state, e, The full measured Pauli set as a function of t, which shows the expected behavior for a CPHASE gate. 
The the y-axes of adjacent elements in the Pauli set are offset by 1. 
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Calibration of RF Sensor Response 

In order to quantitatively interpret sensor values for state to- 
mography, it is important to precisely determine the RF sensor 
response that corresponds to a | S) or a | 7q> state. Because 
the state preparation is imperfect, it is in general difficult to 
accurately measure these values experimentally To provide 
exact calibrations for \S) and |r ), we exploit the fact that our 
sensor is capable of single shot readout. Histograms of sensor 
values for typical measurements yield a double-peaked curve- 
one peak corresponds to \Tq) and one to \S) (Suppl. Fig la). 
In order to calibrate the sensor we first measure T\ at the 
measurement point by preparing a state that is majority |7b) 
(done with a 7i-pulse around the x-axis) and fitting the sensor 
signal to a decaying exponential function of time elapsed before 
measurement (Suppl Fig. la). We note that the measured value 
of T\ is a strong function of the power of the RF excitation 
used to read the conductance of the sensing QD. With prior 
knowledge of T\, we use a procedure similar to that described 
in Barthel et. al (ref. 14) to optimize the measurement time 
given our signal to noise ratios and T\. This process is repeated 
several times per day to check for drift. We recalibrate the 
sensor signals that correspond to \S) and | 7q> for each dataset 
(typically 10 minutes of acquisition time). For each set, we 
prepare a histogram of all observed sensor values. The presence 
of several reference measurements in each dataset guarantees 
that there will be a significant fraction of both \S) and |7b). 
We then fit this double peaked curve to an analytic expression 
corresponding to a weighted sum of two Gaussians with some 
filling in due to T\ decay during measurement (Suppl. Fig. lb, 
purple line) as in ref. 14. From this, we extract the expected 
sensor distributions for \S) and | Tq) (blue and red lines in Suppl. 
Fig. lb, respectively), as well as the fractions of \S) and | Tq) 
present. The centers of the two distributions correspond to the 
sensor signals that will be measured for pure \S) and pure | Tq), 
and using these values we can accurately scale the tomography 
data. We note that this procedure is insensitive to the percent- 
ages of | S) and |r ). In our state tomography only expectation 
values are needed, so the single -shot capability of our readout 
is not necessary beyond this calibration. Nonetheless, we note 
that for the data presented, we measured readout fidelities of 
97% and 98% for the left and right qubits, respectively. 

Readout Crosstalk 

For accurate state tomography it is important that the read- 
out of the two qubits be independent. The two-qubit coupling 
relies on the fact that the \S) and \Tq) in one qubit electro- 
statically gate the other qubit, and reading out the two qubits 
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Supplementary Figure 1 | Singleshot Readout: a, The difference in sensor 
signal between |5>and \ Tq) is fit to a decaying exponential to determine T\, 
which is used in calibration of sensor values, b, The histograms of a mixture of 
|S) and | Tq) states used to calibrate the sensor values. If we choose a threshold 
y thresh t0 distinguish between \S) and |7o) we see a readout fidelity of 97%. 
Purple: fit to noisy distribution including T\ decay from | To> to \S). The 
deduced distribution for \S) (blue) is a Gaussian, while that for \ T$) (red) has 
a tail due to T\ decay. 

simultaneously leads to readout crosstalk; the left- qubit sensor 
value will be different for \S) \ To) and | S) \S), and similarly for the 
right qubit (Suppl. Fig. 2d-e). We avoid this problem by reading 
out the qubits sequentially: while one qubit is read out the other 
qubit is "parked" deep inside of (0,2) where the (0,2) | Tq) is lower 
in energy than the (1,1) | To) so that both the |5)and | To) occupy 
(0,2) (Suppl. Fig. 2a- c). In this way each qubit is read out while 
the other qubit has the same charge distribution for the both 
qubit states. T\ in this region (with the RF excitation off) is large, 
so no measurable degradation of the qubit being stored occurs. 

Calibration of State Tomography 

State tomography involves reading the projection of a qubit 
on to the three Cartesian axes of the Bloch sphere 15 . However, 
the charge sensor next to each qubit allows us to determine only 
the projection on to the z-axis (S-Tq) of the qubit. Therefore, 
in order to perform state tomography, we apply rotations that 
map the x and y axes of the qubit to the z-axis 11 . The y- 
axis component is mapped on to the minus z-axis by a nl2- 
pulse around the x-axis, which is driven by AB Z , and the x-axis 
component is mapped on to the z-axis by adiabatically turning 
on /, which maps the eigenstate | t|}(| |f» into \S) (|7o» (Fig. 
2a). 

Extreme care must be taken when performing two-qubit state 
tomography, as errors in tomography can introduce spurious 
correlations which might inflate the measured degree of en- 
tanglement. However, in S-7o qubits, traditional methods for 
calibration are difficult to implement. The two control axes, / 
and AB Z , are not orthogonal, and precise 90 degree rotations 
are problematic due to the timing resolution available on our 
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Supplementary Figure 2 | Readout Crosstalk: a, A schematic of the energy diagram as a function of e that describes the two qubits and shows the regions of e 
where different operations are carried out. b, A schematic of the signal from the RF charge sensor as a function of e for the two qubits. This signal reflects the 
charge distribution of the two qubit states. For large positive e there is a region where \S) and | Tq) have the same sensor signal (charge distribution), which is the 
foundation of the crosstalk- free readout scheme, c, A schematic of the readout scheme that eliminates crosstalk. First, the left qubit is read while the right qubit is 
"parked" in (0,2), and then the right qubit is read while the left qubit is "parked" in (0,2). d, A two dimensional histogram of the RF sensor responses without (left 
panel) and with (right panel) this crosstalk- free readout scheme, e, Histograms of sensor values without (top) and with (bottom) the crosstalk- free readout. Without 
the crosstalk- free readout the sensor signal of one qubit depends on the state of the other qubit. 



signal generators as well as pulse rise times. Pulse rise times 
in the experimental apparatus prevent instantaneous changes 
in /, which changes both the axis of rotation and the total 
angle rotated. These rise times also cause sudden changes 
in / to become somewhat adiabatic. Furthermore, the effects 
of these pulse rise times depend both on the starting and 
ending point of the an individual pulse. Additionally, the two 
rotations involved in tomography of S-Tq qubits are not easily 
tuned: the accuracies of both the x and y-axis readouts are 
limited by the speed of the waveform generators (~lns in this 
work). Therefore, we perform careful measurements in order 
to determine the three axes on to which the qubit is projected, 
and apply a transformation in order to map these axes onto the 
traditional Cartesian axes. Dephasing during readout rotations 
is represented in this process as readout axes that are longer 
than the radius of Bloch sphere. 

The procedure for determining the the tomographic axes 
makes minimal assumptions. We first calibrate the the S-r 
readout using the singleshot histograms (see above). We then 
assume that any path of the state of the qubit around the Bloch 
sphere in a free-induction-decay experiment should smoothly 
dephase, i.e., there should not be oscillations in the amplitude 
of the Bloch vector. No assumptions are made about the axis 
or frequency of any rotation. If either the length or the angle 
of one of the tomography axes is incorrect, we expect to see 
"ripples" in the length of the measured Bloch vector as a func- 
tion of evolution time (Fig. 3b -c). If we simultaneously consider 
many paths of evolution around the Bloch sphere, errors in 
the angles and lengths of the different readout axes become 




Supplementary Figure 3 | Calibration of state tomography: a, A schematic of 
how state tomography is performed in S-Tq qubits 11 . The S-7o component is 
read out by charge sensing in a region where \ S) occupies (0,2) and | Tq) occupies 
(1,1). The x (y) component is read adiabatically turning on / (rotating by jt/2 
around the x-axis), followed by charge sensing, b, A schematic of the length of 
the Bloch vector for perfect state tomography, c, A schematic of the "ripples" in 
the length of the Bloch vector if state tomography is flawed. The tomographic 
axes are determined by minimizing these ripples. 

distinguishable due to varying phases, periods, and amplitudes 
of the ripples in the lengths of the measured vectors (Suppl. 
Fig. 4a). Therefore, to calibrate our axes, we gather data on 




50 100 150 

Rotation Time (ns) 

IS) 



50 100 
Rotation Time (ns) 

IS) 



150 




o 
-i 



-1 
c i 



-1 







_■ 


d 


Raw Data 


■ 






■ ■_ 


_■ _ 


_■ 


■ 


Rotated Data 


■ ■ 


■ 


■ 


■ 


Target State 



<XI> 
<YI> 
<ZI> 
<IX> 
<IY> 
<IZ> 
<XY> 

<xz> 

<YX> 
<YZ> 

<zx> 

<ZY> 
<XX> 
<YY> 

<zz> 



Pauli Set of Raw Data 

<XI> 
<YI> 
<ZI> 
(IX) 
<IY> 
<IZ> 
<XY> 

<xz> 

<YX> 
<YZ> 
<ZX> 
<ZY> 
<XX> 
<YY> 

<zz> 



Left Qu bit 




400 800 1 

T(ns) 

Pauli Set of Rotated Data 



200 400 



600 800 1000 1200 
T(ns) 



200 400 



600 800 1000 1200 
T(ns) 



Supplementary Figure 5 | Single-qubit rotations: a, The Pauli set for r = 100ns 
as measured is complicated by single qubit rotations, b, Numerically rotating 
each qubit around the S- Tq -axis simplifies presentation and analysis, c, The 
expected state for r = 100ns. d, The single qubit rotation angles for both qubits 
as a function of t are smooth and monotonic functions, e-f The entire Pauli set 
as a function of t for the raw and rotated data equation(l). The the y-axes of 
adjacent elements in the Pauli set are offset by 1. 
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Supplementary Figure 4 | Calibrated vs uncalibrated state tomography: a, 

Data taken to calibrate the tomography shows ripples in the length of the Bloch 
vector if we assume that the tomography projects the quantum state on to 
Cartesian axes (inset), b-c, The paths around the Bloch sphere for the different 
evolutions that are used for tomography calibration. If the tomography is 
assumed to project on to the Cartesian axes there are points that lay outside 
the Bloch sphere, and the pure \S) states are not at the north pole, which is 
indicative of flawed state tomography, d, The ripples in the length of the bloch 
vector are diminished (compared to panel a) if the axes deduced from state 
tomography (inset) are used, e-f, The paths around the Bloch sphere for the 
different evolutions that are used for state tomography. When the correct axes 
are used, all the points lie inside the Bloch sphere and the pure \S) are at the 
north pole. 

many different evolutions around the Bloch sphere by evolving 
from many different starting points at many values of e (Suppl. 
Fig. 4b,c,e,f). We determine the axes on to which we project 
our state by finding the axes that minimize the amplitude of 
the ripples in the length of the Bloch vectors (Suppl. Fig. 4d). 
Based on our measurement procedure, we define the S-T axis 
to lie along the z-axis. We allow the y-axis to lie anywhere 
on the Bloch sphere because a rotation around the x-axis can 
suffer from over/under rotation as well as adiabaticity issues 
with switching / on and off instantly. We constrain the x-axis 
to lie in the x-z-plane because the only expected error is due 
to adiabaticity turning / on and off. The typical tomographic 
axes are shown in Suppl. Fig. 4d, and the signs of the errors are 
consistent with their origins. The variation from calibration to 
calibration is -1% on the axis lengths and angles. 



Determining Single Qubit Rotations 

During the entangling sequence the two qubits rotate very 
rapidly around the S-Tq axis compared to the speed of the 
CPHASE gate {h\m ~ h\m ~ 300MHz, hilm ~ 1MHz). These 
single qubit rotations are not perfectly canceled out by the n- 
pulses in the dynamically decoupled sequence due to pulse 
distortions, consistent with pulse rise time effects at short times 
and capacitive coupling to RC-filtered DC gates at long times. 
Moreover, the angles by which the qubits are rotated change 
as a function of the evolution time t. In order to undo these 
rotations, we perform a least- square fit of the data to the ex- 
pected form of the Pauli set (see equation(l) below), restricting 
the rotation to be around the S-7 axis because 7i,/2 » &B Z . 
These angles are shown in Fig. 3b, and exhibit a smooth, 
monotonic behavior. The angles increase quickly for small r, 
which is consistent with pulse rise time effects, and display 
linear behavior for long r, which is consistent with long time RC 
filtering. For comparison, we plot the entire Pauli set for both 
the rotated and unrotated data in Suppl. Fig. 5 c-d. 



Fit Form for a Dephased Bell State 

In order to fit the experimentally measured Pauli sets we 
calculate p {t) in the presence of fluctuations on J\, and h, and 
average over all these fluctuations 25,26 . For the present case 
where J\2 « /i, h > we neglect the two-qubit dephasing, which 
is smaller than single-qubit dephasings by a factor of j^, « 
300. This yields the following non-zero elements of the Pauli set 



where T 2 j are the single-qubit coherence times of the two 
qubits. From fits of the data to this form we extract for 8e =0: 
T 2> \ =420ns, T 2>2 =510ns, and Ju/2n =0.87MHz. Imperfect state 
preparation, in which the 7r/2-pulse does not leave the state 
of the qubits in the x-y-plane causes mixing between terms, 
which is visible in the form of small amplitude oscillations in 
the {XI), {IX), and {XX) components of the Pauli set (Suppl. 
Fig. 5f). 
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e~ rlT2 ' 1 cos(/i 2 t) 
e~ rlT2 ' 2 COS(/i 2 T) 
e~ r/T ^ sin(/ 12 T) 

e~ T,T2 ' 2 sin(/i 2 T) 
e ~ r/T 2,i e - r /T 2 ,2 
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